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be used to study universal properties of time-reversal invariant disordered systems 
which are either insulators or close to the Anderson localization transition. We derive 
a two-level contribution to the correlation functions of the generic ADRMT and apply 
these results to the critical (multifractal) power law banded ADRMT. Analytical results 
are compared with numerical ones. 
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1. Introduction 

Random matrix tlieories (RMT) have proven to be useful mathematical models to 
study quantum mechanical disordered or chaotic physical systems pp. The Wigner- 
Dyson RMT (WDRMT), which is characterized by independent Gaussian probability 
distributions of the matrix elements with constant variances of the off-diagonal matrix 
entries, can be used to describe the statistics of electrons in small disordered samples in 
the universal (metallic) limit, as the supersymmetric a-model obtained for the WDRMT 
is the same as that of electrons in a metallic dot with infinite conductance [21 E]- The 
ensembles in the WDRMT can be classified by their symmetries, for example, the 
Gaussian orthogonal ensemble (GOE) consists of real symmetric matrices and can be 
used to describe time-reversal invariant systems; the Gaussian unitary ensemble (GUE) 
consists of Hermitian matrices with equal variances of the real and imaginary parts 
of the off-diagonal matrix elements and can be used to describe systems with broken 
time-reversal symmetry [1]. 

Beside the universal metallic regime, physical applications require RMT models 
which can be used in cases of either insulators or critical systems at (or close to) the 
Anderson localization transition [5]. In such unconventional RMTs, the variance of the 
matrix elements depends on their distance to the diagonal. For example, the power law 
banded RMT (PLBRMT) is an unconventional RMT where the off-diagonal variances 
decrease as a power law r"^"' of their distance to the diagonal r outside of a band of 
width B [6]. For a = 1, the PLBRMT exhibits multifractal eigenfunctions and critical 
level statistics, a behaviour different to the Poisson statistics of ideal insulators and 
the Wigner-Dyson statistics of ideal metals [H]. Numerical studies reveal an excellent 
agreement between different correlation functions of the critical PLBRMT and of the 
Anderson model at the transition point [H El E] • The bandwidth controls the fractality 
of the RMT eigenstates: their fractal dimensions are close to the space dimension, d, if 
the bandwidth is large and are much smaller than d in the case of the almost diagonal 
critical PLBRMT with small bandwidth. The latter regime of strong multifractality is 
relevant for the Anderson model in high dimensions. 

The methods which can be used for an analytical study of the PLBRMT 
depend on the bandwidth. For large bandwidths, a field theoretical technique of the 
supersymmetric a-model can be applied [6j. In the opposite case of the almost diagonal 
PLBRMT with small bandwidth, a complementary method, the supersymmetric virial 
expansion (VE), can be used. The VE is a regular expansion of the correlation 
functions in a number of interacting energy levels. In the supersymmetric version, 
it corresponds to the formal expansion in a number of independent supermatrices with 
broken supersymmetry. It is the supersymmetric analogy to the linked cluster expansion 
used to calculate thermodynamic quantities for classical imperfect gases [lOj. 

The VE has been initially developed on the formal basis of the Trotter formula 
and a combinatorial analysis for a spectral form-factor of the generic case of almost 
diagonal random matrices (ADRMT) from the GUE symmetry class [11]. It has then 
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been extended to calculate the density-of-states and the level compressibility [121 IE]- 
The supersymmetric (SuSy) VE has been developed later for unitary ADRMTs |14j . 
It is free of the complicated combinatorics and allows one to investigate not only the 
spectral statistics but also the statistics of the eigenfunctions. Recently, the SuSy VE 
has been used to study the dynamical scaling in multifractal systems [15]. 

To use the SuSy VE in the context of physical systems with time-reversal invariance 
it is necessary to modify it such that it is applicable to orthogonal ADRMTs. Many 
experimentally relevant disordered quantum systems, e.g. cold quantum gases in a 
disordered optical lattice, belong to this class of systems [IHl CZ] . This paper fills this 
gap in the theory of the ADRMTs, namely, we extend the general formalism of the 
supersymmetric VE to orthogonal ADRMTs. 

We organise the paper as follows: basic definitions of the almost diagonal RMT 
and of the two point correlation functions are given in Section 2; the supersymmetric 
representation of the correlation functions and of the supersymmetric VE are explained 
in Sections 3 and 4 respectively; a parametrization of supermatrices is given in Section 5; 
generic results for the leading terms of the VE are calculated in Section 6 and are applied 
to study the local density-of-states (LDOS) of the almost diagonal critical PLBRMT in 
Section 7. To the best of our knowledge, the LDOS-LDOS correlation function of the 
critical orthogonal PLBRMT was never calculated before. We compare the analytical 
results with the numerical ones and conclude the paper with a brief discussion. 



2. Main Definitions 



We consider a Gaussian real symmetric ADRMT defined by the following statistics of 
matrix entries: 

{H,,) = 0, {Hl) = l, (ijy = i?2F(K-j|)=6,,«l, (1) 

where F{x) is a real decaying function and the small parameter B <^ 1 reflects that 
the ensemble is almost diagonal. The matrices of the ensemble are of size N x N with 
N ^ 1. We introduce the retarded and advanced Green's functions by 

G^^^'iE)^ I ^-+0; (2) 

E — H ± IT] 

with H being an arbitrary matrix of the ensemble and we define the two point Green's 
function at the band center £" = as 

- (I) Gi (-|) , (3) 

using the matrix elements of the Green's functions Gnn^ = (nlG^^'^ln) ( \n) with 
n = 1 . . . is a canonical basis vector in the vector space of the RMT ) . 

We will consider disorder averaged two point-correlation functions, namely, the 
spectral correlation function and the diagonal local-density-of-states (LDOS) correlation 
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function, which are defined by 

N 

R,i^)^A^Y.{{s[E, + ^)s{E,-^))), (4) 

and 

C,{u) {{\^^iP)\" l^MfS [e, + I) 8 {e, - 1))^ . (5) 

p=i ij=i 

Here A is the mean level spacing of the RMT; Eij are the eigenvalues of the random 
matrices; {{ab)) = {ah) — {a){b) and (...) denotes the average over the ensemble of 
random matrices. For orthogonal almost diagonal RMT, the mean level spacing is 
known to be [TTj: 



N N 



(6) 



The spectral and the LDOS correlation functions can be expressed by the two point 
Green's function: 

A2 ^ 

^2M = — $^Re[((^p,H))], 



27r2 



A2 

271' 



p,q=l 



J]Re[((^,,M))] 



(7) 
(8) 



p=i 



Thus, one needs {{Qpqioj))) to study the correlation functions 



3. Supersymmetric Representation of the Two Point Green's Function 

The ensemble averaged two point Green's function can be obtained with the help of a 
supersymmetric field theory [H]. Let us introduce 2N super vectors of the form 

R/A 

(9) 



R/A 



\ 



sry 

R/A\* 

X. ) I 



t 



where i = 1, . . . ,N, Sf^^ are commuting and xf'^ anti-commuting Grassmannian 
variables and the indices R/A stand for the retarded/advanced sectors. We use the 
outer product of these super vectors to define the supermatrices 



R/Ay 



(10) 



R/A 



Q^ 



1 

2 



11) 
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These definitions allow one to write the ensemble averaged two point Green's function 
as 

{g„{u;)) = I D{Q}n,A, |^ne^°[«'lj e^^[«-«^l j . (12) 

Here / D{Q} = YliLi D{Qi} is the measure of integration over the supermatrices Qi and 
TZp, Aq are the symmetry-breaking factors between the commuting and ant i- commuting 
variables in the retarded and advanced sectors, respectively. They can be chosen as 
follows 

n, = x^{xi)\ A = x^W)*. (13) 

The action of the field integral is split in two parts: 5*0 depends on a single supermatrix 

S,m (Str[g,])^ + ^i^Str[AQ,], A = |^ J j , (14) 

Sp depends on the product of two supermatrices and is proportional to the off-diagonal 
variances 

5'p[Qi,(5j] = -6ijStr[Qi(5j]. (15) 



4. The Virial Expansion 

The two point Green's function ([T^ of the ADRMT can be approximately calculated 
using the virial expansion [TH [121 [ll]- The basic concept of the VE is similar 
to the idea of resonant level interactions which has been used in a semi-empirical 
renormalization group approach [19]: a small probability for energy levels of the 
ADRMT to interact in the energy space allows one to develop a regular perturbative 
expansion in the number of interacting levels. Each term in the VE reflects a 
contribution of a certain number of simultaneously interacting levels. In contrast to 
the semi-empirical renormalization group approach, the VE is a regular perturbative 
expansion, thus allowing to include an arbitrary number of resonant level interactions 
and a controlled estimate of the contribution of the unused resonances. 

The mathematical expressions of the virial expansion can be obtained by rewriting 
the term containing the off-diagonal variances in the superintegral ( |T2l) 

N oo 

JJ g5p[Q,:,Qj] ^ yi3 ^ ^ yM ^ ^^^^ 
i^j m=2 

N 

= 1 , V(2) = (e2^-[Q-«^] - 1) ; (17) 
i>j=i 

and so on (see details in [Ej). 

The VE is similar to the cluster expansion in the classical gas theory [lOj, however, 
the SuSy representation of the resulting terms ensures that, in contrast to the cluster 
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expansion in the classical gas theory, no further reordering of the terms is necessary. 
Thus, corresponds to the contributions of the diagonal part of the RMT, V*-^-* 
reflects the contributions of two-level interactions and so on. Using Eqs. ( |T6lfT71) . we 
can formally expand the averaged two point Green's function: 

After integrating out all supermatrices Qj which neither included in the SuSy breaking 
factor (j 7^ p,q) nor are linked to this factor through V*-™'-' functions [Ej, the first two 
terms of the VE read 

(g;^^) = J D{gp}D{gJ7^pAe''»[«-l+^»[Q'l, (is) 

i^pl) = J ^{Qp}^{QJ^pA (e^«[«^l+^«[«^]) {e'^^lQ^'^"^ - 1) , (19) 

N 

5. Parametrization of the Supermatrices 



(20) 



A convenient parametrization of the Q-matrices for the SuSy VE has been first suggested 
in ^20j. For each supermatrix, one introduces 2 positive variables Xf^^ G [0, oo], 2 angles 

(pf^^ G [0, 27i] and 4 anti-commuting Grassmann variables rjf^"^, (^rif^"^^ {i = 1, . . . N) 
in the following way: 

r,R/A .R/A f. 1 / R/A\* R/a\ j<A«M R/A .R/A R/A i<h^/^ x 

Si' =\' ( 1 - 2 (^^i j ) e ' ^^ =\ V^ . (21) 

The integration measure in this parametrization is given by 

r. r^. A dX^dcPf {dvn*dv^ dXfd^f (dyf)* dyf 

^^^^-n TTAf irXf ' ^^^^ 

1=1 ' ' 

and the single matrix part of the action reads 

Som = -\ ((Aff - (Xffy + {(Xff + (Xff) . (23) 

Using rotations in superspace, we can find a convenient representation of the two matrix 



part of the action in the new variables (see Appendix A for details). 



-AfA/cos(^j)) (l-^K) 



A\* „A 



) • (24) 



, R/A R/A R/A 1 

where a^j = Vi ~ Vj ^^'^ 

»r - *r - - 5 ( (of") * or - ' of") . (25) 
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As a direct consequence of the underlying symmetries of the ensembles, the 
supermatrices defined in Eq. ffTTj) are twice the size of the supermatrices which have been 
used for the SuSy VE of the unitary ADRMT ( cf. Eq. (A.1),(A.2) [H] ). Parametrizing 
the supermatrices as defined in Eq. (1211) . this yields a different amount of effective 
variables necessary for the integral representation of the two point Green's function; the 
two matrix part of the action contains one effective angle in the case of the GUE ( cf. 
Eq. (A. 10) [14J), whereas for the GOE two effective angles {9jj and 9^ in ([211)) appear. 
For the further calculations, we will use other commuting variables which allow one to 
simplify the single matrix part of the action: 

R. ^ {X^Y ~ {Xt)\ 5.^(Af)%(Af)^ (26) 
where —Si < Ri < Si, Si G [0, oo]. The integration measure changes to 

n r^i _ -rr dsm^d<p^d<j>f [d^^Yd^f {dyfYdyf 

6. Generic Results for the First Two Terms of the Virial Expansion 



Integration of the superintegral in (1181) and calculation of the diagonal contribution of 
the product of single averaged Green's functions (G^ (^^/2) ) {—uj/2) ) yields: 



((CM»=0, ((6°(^)))|^^^^«-H^ = 2x/V*W + H|^, (28) 
where s = a;/A. 

The superintegrals (fT9l) . (120|) can be calculated using a saddle-point approach 
(cf. the Section 4.1 in [E]). Rescaling all Si and Ri variables with Si = ^/h^Si, 
Ri = \/b^Ri, where a e {n, q}, we obtain a large negative factor of order 0(1/ B"^) in 
front of R^ in the single matrix part of the action 

So[Q^] = -7^R■+^p7=S^, (29) 



2&pa 2\fbp(^ 
This allows to employ a saddle point approach in the Ri. 

'"-J-sJ "^{Sl-Rl) {Sl-Rl) K 







02 C2 V / \ ^ P J \ Q 



where / [Rp^a, Sp^a) = [e'^^p^^P'^"^ — l) and Rt, St are the typical values where 
the integrals over R and S converge. Based on rough estimates, one obtains two 
requirements for the validity of the saddle point approximation [T3j: 

< 5 < 1 . (30) 

However, the analytical results derived from the saddle point approximation and the 
results of the numerical diagonalization are in a very good agreement even in the regime 
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cij ~ 1; see an example of the LD OS-correlation function C{uj), Eq. ( H2l) : figure [1] 
below, and figure 3 in the paper [7] in the cases of the orthogonal and the unitary 
critical ADRMT, respectively. Corrections beyond the saddle point approximation will 
be studied in details elsewhere ^2T\ . 

After using the saddle point approximation, we can integrate out the angles, the 
Grassmann variables, and the remaining real variables. This yields 



k=l 



i^sm^f±i-^r{^' ''^-J^^-'\ (32, 

These power series can be summed up by using a Fourier transform of the physically 
relevant real part of the two point Green's function: 

Re [g,M = —j due-^-' {g„ (u) + g;^ (^)) . (33) 

After applying the Fourier transform to Eqs. flHTj) . (1321) . we find the two point Green's 
function in the time domain: 

27r2 
A 

Re [{{gsm = -^E^-i^i^^'^"*"'"'^o M . (35) 

The energy representation of the two point Green's function can be obtained by the 
inverse Fourier transform. If 7^ we find: 



Re 



Gplm = ^&P.|t|e-^-*^-IV, (6,/) , (34) 



Re 



= (I)' 1^1^""^ (3^0 (-^) +/i (-^)) - vr, (36) 



N 



3 



Re mjm] = E (i)' l-l-"'^ (^0 (-^) - h (-^)) ; (37) 

where u = u / {A^Jh^)^ and a = g, n in Eqs. ( l36l) and ( 137|) . respectively. These are the 
generic results for the two matrix approximation of the two point Green's function for 
ADRMTs of the GOE symmetry class. 

7. Application to the Critical PLBRMT 

The critical PLBRMT is defined as follows: 

(l/,> = 0, W5> = 1, = (38) 

The ensemble fl38|) is named "critical" because it shows multifractal behaviour of the 
eigenstates and critical level statistics at arbitrary B 
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7.1. The Spectral and the LDOS Correlation Functions 

Inserting the variances ( !38l) into Eqs. (El), ( !36l) . ( 1371) . we obtain the spectral and 
the LDOS correlation functions for the critical ADRMT of the orthogonal symmetry: 

1 ^ 

^2M = ]^ E {V2^mP - n\e-^'\^-\' lo {^'\p - n\') - l) , (39) 

here lj = uj/{\^B). Let us introduce a shifted (positive) spectral correlator: 

1 ^ 

/?+(.;) ^— J2 V2^mp-n\e-'^'\^~-\'lo{Co'\p-n\') , (41) 
and scale the LDOS correlation function: 

^<-) - Im- (^^' 

This scaling eliminates the level repulsion effect at o; < BA, therefore, C(u;) is 
convenient to study the correlations of the eigenfunctions taken at the same space point 
and with a given energy difference. The unitary case has been explored in The 
qualitative behaviour of C{uj) for the orthogonal case is similar to the unitary one, see 
figure [H There are three regions which reflect different physical phenomena: 

1) uj > Eq oc B: the energy difference is larger than the hopping band width of the 
ADRMT which results in anticorrelations of the eigenstates and in a generic power law 
behaviour C^u) oc cu"^ (see [7] for details). 

2) < uj < BA: the energy separation is very small so that the behaviour of 
6*2 and R2 is governed by the level repulsion in the energy space, therefore, the ratio 
C = C2/R2 saturates to a constant. Note the analytically obtained constant is a bit 
larger than the numerically found value because of the reasons explained below. 

3) BA < UJ < Eq-. this is the most interesting region where strong spatial 
correlations of the fractal eigenstates must result in a scaling C{lj) oc uj'^ , = d2/d — 1, 
d2 and d are the second fractal dimension and the space dimension, respectively [221 ES]- 
This property is especially nontrivial in the case of strong multifractality d2 d 
{d = 1 for RMT) where the fractal eigenstates are sparse, nevertheless, they are strongly 
correlated in space. The fractal dimension ^2 depends on the symmetry class and on 
the parameter i?, for example: d2 = ^/2B for the GOE and d2 = 7r5/v^ for the GUE 
critical ADRMT ( cf. [24J ). Thus, the critical ADRMT corresponds to the case of the 
strong multifractality. The parameter Eo can be referred to as "the upper cut-off of 
the multifractality". It can be shown from Eq. fl40l) that 

Eo ~ ^NAd2 . (43) 
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Figure 1. (coloured online) LDOS correlation function for the critical ADRMT with 
B = 0.06. Main panel: Comparison of the analytical results for the cases of GOE 
and GUE , N — 2000. Black solid lines indicate a reference slope predicted by the 
two matrix approximation of the VE, Eq. ([l5)) Inset: Comparison of the analytical 
and numerical results for the case of GOE, N = 10000. Note that the difference 
between the numerically and the analytically found slope in the intermediate regime 
BA CO <^ Eq is due to the accuracy of the two matrix approximation, see detailed 
discussion in the main text. 



Note that the same result holds true for the GUE ensemble: substituting the unitary 
values of A and ^2 into Eg. (143!) we find the unitary value of Eq [7]: 



Er 



GUE 



B> Er 



GOE 



(44) 



If d2 is small, the leading term of the VE cannot distinguish the exponents fi = d2 — l 
and /i ^ — 1. One has to take into account the interaction of a larger number of energy 
levels to find the correct exponent \T^. That is why the analytical result (1401) shows the 
scaling 

771 

NC(BA < < Eo) ^ — (45) 
see figure [H and the analytically calculated constant for C{uj < BA) is a bit larger as 



compared with the numerically obtained value. As Eq 



> Er 



GUE 



the GUE curve 



GOE 
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in figure [His shifted upwards compared to the GOE curve. 

8. Conclusions 

In this paper we have developed the supersymmetric virial expansion for the two point 
correlation functions of the almost diagonal Gaussian Random Matrix ensembles of the 
orthogonal symmetry class. We have derived the generic results for the two matrix 
approximation of the two point Green's function in the time (Eqs. flH^ . flH3|) ) and 
energy (Eqs. ( l36i) . ( 1371) ) representations. This contribution to the Green's function 
results from an interaction of two energy levels in the energy space and can be used 
to study statistical properties of time-reversal invariant disordered systems which are 
either insulators or close to the point of the Anderson localization transition. 

To demonstrate how the method works, we have applied the generic results of the 
two matrix approximation to the Gaussian orthogonal critical power law banded RMT 
with small bandwidth, Eq. (IHHj) . We have obtained for this ensemble the spectral 
correlation function (Eq. (1391) ) and the local-density-of-states correlation function (Eqs. 
fHUjl . fH^ ). the latter has to the best of our knowledge not been derived before. A 
comparison of the analytically obtained function C{u!) with the results of the direct 
numerical diagonalization shows an excellent agreement at u > B, see the inset in 
figure [U In the region BA < uj < B the two matrix approximation yields qualitatively 
correct results: we find a scaling C{uj) oc uj^^ in this parametrically large energy 
window, however, the exponent is different from the exponent of the expected critical 
scaling C{uj) oc uo'^'^~^ This is because the second fractal dimension is small in the 
case of the ADRMT, ~ -B <^ 1 [24j and the leading term of the VE is unable to 
distinguish exponents 1 — d2 and 1 [7j. One has to take into account the interaction 
of a larger number of the energy levels to find the correct exponent [15]. Nevertheless, 
a coefficient in the scaling is proportional to d2 : C{uj) ^ Ad2/{2Lj), see Eq. f HHI^^j) . 
Thus, even the results of the two matrix approximation refiect multifractality of the 
critical eigenstates [TS]. A Fourier transform of C{uj) yields the return probability for 
an initially localized wave packet. Therefore, the results of Section 7 might be relevant 
to study fractal properties of cold atoms in a disordered optical potential which does 
not break the time reversal symmetry [ini [T7] . 

A generalization of the current results for the three matrix approximation is 
straightforward (cf. [E]). In the future, we plan to use a SuSy field theory, which is 
the starting point of the of the perturbative VE (Eq. ( fT2l) ). to obtain nonperturbative 
results accounting for the interaction of all energy levels at least with an accuracy of 
the first loop of the Renormalization Group procedure. 
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Appendix A. Simplification of the Two Matrix Part of the Action 



U 



\ 



\ 



RA 



R/A 



V 



R/A 



(1+1 (^i?M)*^iiM) 



cr is a symmetric orthogonal matrix 



a = 



a' 
a' 



a 



1 



^ 



1 


1 








1 / 



V2 

\ 

and D is a matrix which is diagonal in each sector 

/ 1 \ 



D 



XrXr —XrXa 
XaXr —XaXa 









(A.i; 



The supermatrices ffTTl) written in the parametrization can be diagonalized in each 
sector by the following transformation: 

D = a'^U^QUa-^ , Q = UaDaUl 

Here U and W are unitary block diagonal matrices: 

/ Urc''^'' 

U*Re-"f''' 

Uac'^^ 

UXe-''^'' J 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



In this representation R labels the retarded and A the advanced sectors of the 
supermatrix. Each of the 4x4 submatrices of the sectors in Eqs. (lA.4l) .( lAr5l) contains 
bosonic and fermionic sectors in a pattern determined by the outer product of the super 
vectors defined in Eq. ( JTOl) . The supertrace in the two matrix part of the action can be 
simplified by using the invariance of the supertrace under cyclic permutations 



Str [Q,Q,] = Str 



Str 



aUpiaDiaUjUjaDj 



UiaDiaUjU^aDjaUj 
and the following property of 2 x 2 supermatrices (cf. [25 

This allows one to obtain Eq. f l24|) . 



(A.6) 
(A.7) 
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